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Tftte  ftpovt  focuses  cm  s  USGS  provkied  diMa  set  from  refractkm  siirve3ni 
wMdn  ttie  State  ot  Makie.  These  ckta  contain  excdtoit  of  siKVt 

period  Eayfcigh  smves  (Rg)  generated  hy  the  explosions.  The  iii^at  inters 
was  to  ose  tfiese  data  to  characterize  shear  wave  vdodty  and  Q  ki  the  ixf^per 
1>2  km  in  Ms  envaronsacnt  of  fokied  palecxraic  secMients.  However  die  chita 
set  csdAitted  imeresting  wave  profx^atkm  charactertetks.  Tim  principle 
rcsidls  of  die  study  are  as  fbSow: 

Rnt,  s%Rtficant  lateral  duuiges  in  die  moveout  of  the  Rg,  S  and  P  wave 
were  observed.  Anaijrsis  techniques  were  devekqied  to  successfufly  invert  for 
a  lateraBy  varying  earth  structure  and  to  mexid  the  Rg  wavefmins,  both  in 
terms  of  absedute  amplitude  and  phase. 

Second,  the  shear  wave  vdodty  modd  derived  from  the  surface  waves, 
was  used  in  a  bMnd  test  to  jxredkt  the  S-wave  arrival  times  for  two  pre^es 
havii^;  significant  S  waves.  An  excdlent  fit  was  obtained  without  havir^ 
used  the  S-wave  data  in  the  invarskm  fm  the  earth  modd.  This  indicated 
die  abffity  oi  the  1-4  Hz  surface  wave  data  to  resedve  shear  wave  velocity 
structure  to  depths  of  2  km  in  this  environment. 

Hnalty,  sinoe  absedute  instrument  gains  were  known,  forward  modding 
was  used  to  test  the  laterally  varying  vdkxaty  and  Q  modd,  in  order  to  com¬ 
pare  the  observed  and  estimate  amplitude  decay  in  difrerent  frequency 
bands.  This  permitted  an  estimation  of  the  isotropic  moment  fexr  them  one 
ton  exfdoaknis.  By  comparing  estimates  made  in  different  portions  of  the  1-5 
Hz  fre^ency  baiid,  we  conciude  that  a  step  source  time  function  in  moment 
is  adequate  to  mcxld  the  surfsce  wave  obs^atkms  fmr  these  cme  ton  shots. 
The  estimated  isotropic  mennents  are  similar  to  other  chemical  yield-moment 
observations.  However,  larger  industrial  chemical  Masts,  e.g.,  100  ton  or 
greater,  would  be  expected  to  show  a  frequency  dependence  in  the  isotropic 
moment  estimate  because  of  the  distributed  source  dimensions. 

These  report  presents  the  data  sets  and  condusions  in  four  sections. 


SURFACE  WAVE  INVERSION  FOR  SHEAR  WAVE  VELOCITY 


R.  B.  HERRMANN  and  G.  I.  AL-EQABI 

Department  of  Earth  and  Atmospheric  Sciences 
Saint  Louis  University 
3507  Ladede  Avenue 
St.  Louis,  Missouri  63103 

ABSTRACT 

Land  based  seismic  exploration  has  typically  viewed  the  surface  wave  or 
ground  roll  as  undesirable  noise  while  earthquake  seismology  has  been  able 
to  use  the  surface  wave  to  invert  for  lateral  variations  in  crust-upper  mantle 
shear-wave  velodty.  Although  stirface-wave  inversion  techniques  for  shear- 
wave  velodty  and  shear-wave  Q  are  easy  to  apply  once  the  phase  velodty, 
group  velodty  and  andastic  attenuation  (topersion  are  known,  the  manner  in 
whidt  these  parameters  are  obtained  elfidently  from  large  exploration  data 
sets  has  not  been  fuUy  explored.  This  paper  will  address  the  processing  steps 
required  and  apply  them  to  surface-waves  recorded  in  a  refraction  experiment 
in  the  State  of  ^ine  in  the  U.  S. 

INTRODUCTION 

In  earthquake  seismology,  surface  waves  have  long  been  used  to  infer  shear 
wave  vdodty  modds  of  the  crust  and  upper  mantle  of  the  earth  (Ewing  et  al., 
1957).  Andysis  techniques  have  been  developed  to  invert  for  ^ear  vdodty 
as  a  function  erf  depth,  and  recendy  progress  has  been  made  in  using  these 
data  to  invert  fox  laterd  variations  in  earth  structure  (Nolet,  1987;  Nolet, 
1990).  On  the  other  hand,  less  attention  has  been  paid  to  surface  waves  in 
seismic  exploration,  where  the  surface  wave  is  viev^  as  undesirable  noise 
that  affects  the  imaging  of  subsurface  features. 

The  objective  of  the  paper  is  to  demonstrate  the  use  of  surface  wave  data 
acquired  as  part  of  a  rdraction  survey  for  the  determination  of  shallow  shear 
velocity  and  Q.  Even  though  surface-wave  data  are  available  in  the  distance 
range  of  0  -  40  km,  which  is  much  larger  than  in  typied  seismic  exploration 
data  spreads,  the  need  to  process  multichannd  ^ta  is  a  shared  feature. 
Andysis  techniques  are  revi^ved,  and  are  actud  data  sets  are  processed. 

DATA  PROeXSSING  FOR  DISPERSION 

Given  multiduuuid  data,  disperdon  information  can  be  obtained  by  sin¬ 
gle  duuind  as  wdl  as  multichannd  im>cesdng.  Single  channd  processing 
can  ]^rovide  group  vdodty  dispersimt,  while  multichannd  processing  is 
re<paired  to  obtain  phase  vdodties  without  making  assumptions  about  source 
phase  or  the  number  wavdengths  between  the  source  and  dngle  reedver. 
The  techniques  of  multipie  filter  analyds  and  phase  vdodty  stadefog  are  used 
to  define  the  disperdon  parametars. 


Fgter  Amfyas 

Multiple  titer  tfta^reit  was  introduced  by  Dziewonski  et  d.  (1969)  for  the 
detannim^ion  oi  group  vdodty  from  a  dispersed  waveform.  Hernnann 
(197^  stndted  die  use  oi  imdt^tie  filter  ana^^  to  estimate  q>ectfal  am]^> 
tudes  oi  die  wurious  modes.  Lrt  the  Fourte  translmrm  of  disposed  j»rO' 
pagating  wave  at  ctislanoe  r  be  Define  a  Gaussian  filter  as 


ttf  ■  «(j(ir/o)*^ 

The  result  oi  applying  this  filter  to  a  dispersed  signal  is  a  filtered  time  signal 
vdiose  envdope  readies  a  maximum  at  the  group  vdodty  arrival  times.  The 
envde^  of  dw  ffltered  time  Itistory  for  a  sin^  mode  is  approximatdy 

Cmnputationally^  the  time  series  is  filtered  using  (1),  a  Hilbert  transform  of 
the  filtered  time  series  is  taken,  and  the  two  are  combined  to  form  the 
envdope.  The  envdope  is  searched  for  the  rdative  maxima,  from  which  the 
group  vdodty  is  computed  firom  the  time  of  the  peak  amplitude  and  the 
specbral  an^tiitude  is  estimated  from  (2)  by  the  rdation 

For  ease  of  analyds  each  two  plots  are  made  fcMr  each  trace  analyzed.  The 
first  presents  the  qiectral  amplitude  vosus  frequency,  annotating  each  with  a 
S3nnM  rdaled  to  the  value  of  the  rdative  amjditude.  The  secemd  plot  pro¬ 
vides  contours  of  die  envdope  as  a  function  oi  group  vdodty  and  frequeiuy, 
with  the  group  vdodties  the  envdope  peaks  indicated  by  the  same  symbol 
as  die  fi^  plot.  Examination  (rf  both  plots  assists  in  the  identification  of 
modes. 

Fear  a  given  shot,  individual  traces  are  processed  using  multiple  filter 
analysis,  and  a  taUe  containing  the  frequency,  group  vdodty  and  spectral 
amplitude  is  produced.  Assuming  a  laterally  homogeneous  waveguide,  the 
individud  group  vdodty  estimates  of  a  given  mode  will  be  identied,  and  the 
decay  of  spectrd  amplitude  with  distance,  r,  wiU  have  a  functiond  form 

A(«,r)  ■  r~^.  (3) 

An  interactive  program,  idiich  wmrks  on  non-graphic  computer  displays,  is 
used  fird  to  identity  the  modes  on  the  basis  of  a  group  velodty  histc^am, 
and  tl^  to  define  the  andastic  attenuation  coeffident,  y,  through  a  least 
squares  fit  of  (Observed  spectrd  amplitude  of  this  mode  according  to  the 
oumM  (3).  Provisions  are  made  for  editing  extreme  outliers.  The  result  of  this 
interactive  program  is  a  tabulation  of  group  vdodties  and  andastic  attenua- 
tkm  values  to^rther  with  their  standard  errors,  for  use  with  an  inversion  pro¬ 
gram. 


Group  Velocity  Stacking 

The  above  technique  for  estimating  the  average  group  velocity  and  ane> 
lastic  attenuation  over  the  profile  requires  running  multiple  filter  analysis  for 
each  trace,  consolidation  of  output  into  a  large  file,  and  then  interactively 
editing  the  data  set. 

Barker  (1988)  introduced  the  idea  of  determining  the  group  vdodty  by 
staddng  the  envelopes  of  different  traces  for  a  given  frequency.  This  is  possi¬ 
ble  if  the  envelopes  are  transf(»med  from  a  function  of  time  to  a  function  of 
group  velocity.  The  group  velocity  of  the  spread  is  obtained  by  searching  for 
Ae  maximum  of  the  stadc.  If  a  program  automatically  selects  the  maximum, 
it  can  then  return  to  the  envelope  of  each  trace  to  search  for  the  maximum  in 
the  neighborhood  of  the  velocities  from  the  automatic  selection  so  that  error 
estimates  of  the  group  velocity  can  be  made.  In  addition,  spectral  amplitudes 
can  be  used  to  estimate  anelastic  attenuation  coefficients. 

Phase  Velocity  Stacking 

M<Mechan  and  Yedlin  (1981)  introduced  a  techiuque  for  the  determina¬ 
tion  of  phase  velocity  by  pi^orming  a  p  -t  stack  followed  by  a  transforma¬ 
tion  into  the  p  -u  domain.  Mokhtar  et  al  (1988)  described  a  procedure  for 
accomplishing  this  by  first  Fourier  transforming  the  time  series,  and  then  per¬ 
forming  the  stacking  operation  in  the  frequency  domain.  Let  the  Fourier 
spectrum  of  a  seismic  signal  of  the  i'th  station  at  distance  be 

A{ta,ri)e^^**^,  (4) 

where  the  phase  t«rm  i|f(a>)  includes  a  distant  dependent  component  for  pro¬ 
pagation.  The  p-u  stack  of  N  traces  is  defined  by  the  relation 

l/(P,«)-  2C(w)-‘i4(a»,r<)e^*'e^'^'.  (5) 

i>l 

where 

C(w)  ■  A(io,ri)e^'^. 

Division  by  C(<u)  is  a  simple  artifice  to  normalize  the  spectrum  and  remove 
the  source  phase  if  the  signal  is  that  of  a  single  mode.  The  observed  spectnun 
is  assumed  to  be  the  superposition  of  M  sui^ce-wave  modes  such  that 

where  5|(«i>,rj)  is  the  amfditude  spectrum  of  the  it'th  mode  and  the  phase  is 
separated  into  distance-independent,  ^(a((o),  and  distance-dependent, 
(0poi((a)rj  components,  pof^  is  the  inverse  of  ffie  phase  velocity  of  the  I^'th 
mode  at  frequency  m. 

If  the  sigiud  consists  primarily  oi  a  sin^e  noise  free  surfiice-wave  mode  at  a 
giv«i  frequency,  ffien  the  quantity  U(p,<a)  will  have  a  maximum  when 
p  «  p^.  Alternativdy  it  is  convenient  to  sear^  for  the  maximum  of  the  func- 
tton 

I  lf(p,»)  I  (7) 
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to  find  the  di^Mersion  curves  since  it  b  a  real  function.  Since  there  are  N  dis- 
IWAfiiir  die  maadmum  i^Uue  of  dte  quantity  i  U(p,e>)  I  should  be  equal  to  N. 

tMs  ii^due  is  less  than  N,  then  we  may  attribute  this  to  an  enm  in  dte  ray 
pimmseter  between  tiKe  stations  and  thus  estimate  the  error  in  the 
ffnfrsipendin^  phase  vdodty.  If  we  assume  that  in  the  neighborhood  of  a 
Stiidk  maximum  die  quantity  Ap  -  p  has  a  normal  distribution  with  zero 
mean  wad  a  vmianoe  dien  the  expected  value  of  any  twm  in  (4)  is 

(8) 

whme  we  used  the  definition  a>  «  Ittf. 

Mnoe  eadi  term  in  (7)  is  always  positive,  the  expected  value  of  the  stadc  (5) 
of  a  sin^  mode  is  just 


Qven  the  stadc  value  (7),  a  Newton>Raphson  technique  is  used  to  find  the 
value  of  <r  from  (9)  that  corresponds  to  this  value.  The  error  in  phase  vdodty 
is  obtained  using  the  definition  c  «l/p  and  the  rdation 


This  rdation  was  tested  by  numerically  modeling  a  staddng  operation  in 
which  the  ray  parameter  errmr  had  the  assumed  disMbution. 

In  the  more  realistic  case  of  multimode  surface  waves,  (7)  will  not  yield  a 
maxinium  independent  of  the  amplitude  spectrum  of  the  other  modes.  Thus 
the  stadc  value  at  a  given  frequency  will  typically  be  largest  for  one  mode 
and  smatter  lor  others.  The  simplistic  error  analysis  will  yield  larger,  perhaps 
unumrranted,  errors  for  the  other  modes.  This  is  an  inherent  problem  with 
this  tedmique,  that  can  only  be  resolved  if  a  phase  matched  filter  tedinique  is 
first  apffiied  to  each  input  ipectra  to  isolate  a  single  mode  befcMre  the  phase 
velodty  stadc  is  perfcnni^. 

Evaluation  of  Technujues 

All  these  techniques  suffer  from  the  need  for  manual  intervention  to 
cmractly  define  the  dispersion  data,  especially  identificaticm  of  the  mode. 
The  on^  real  identification  of  the  mode  numb^  is  to  count  the  zero  crossings 
el  the  aur:bce>wave  eigenfunctions,  with  the  fundamental  mode  Love  and 
vertical  coa^ponent  Rayldgh  e^enhinctions  having  zero  zero-crossings  with 
depth,  fiigenfunction-depth  data  are  not  available  with  surface  sensors.  An 
ahemative  technique  uses  the  fact  that  the  fundamental  mode  has  the  lowest 
phwWL.  amd  usualhc  poup,  vdodties  at  any  frequency.  For  swface  sotuces, 
die  fiuwiaiatrdai  mow  is  often  the  one  with  largest  ampfitudes. 

The  fiiel  that  the  ^>ectral  amplitude  level  of  the  surtece  wave  anodes  varies 
wMik  fircqaency,  the  source  exdtatkm  fcxrces,  and  with  distame  doe  to  aadbs- 
tlc  attenwalinBi,  aacana  that  not  necessarily  the  smxw  mode  is  deaaiwawt  a4  a 
dews  fiequwncy.  Because  of  this,  the  plu»e  velodty  stacking  and  ancdd^e 
Wtet  tedadqucs  wffl  find  it  difficult  to  follow  the  dispersion  curves  of  a  low 
Wiwlitttdt  mode  in  the  presence  of  a  larger  amplitude. 

Tm  gproup  vdodty  stacking  technique  %vorks  reasonably  well  for  estimating 
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the  group  velocity  value,  even  though  the  determination  of  the  average 
^oup  velocity  should  be  weighted  as  a  function  of  distance,  since  the  group 
velocity  resolution  is  less  at  shorter  distances.  The  automatic  determination  of 
the  anelastic  attenuation  coefficient,  7,  also  requires  such  a  weighting.  In 
addition,  ^ce  this  depends  upon  spectral  amplitude  estimates,  the  interfer¬ 
ence  of  modes  at  short  distances  will  give  erroneously  spectral  amplitudes, 
and  thus  overestimate  the  y  values.  In  such  a  case,  it  is  better  to  run  the  mul¬ 
tiple  filter  analysis  on  each  individual  trace,  and  then  interactively  determine 
the  values. 


SURFACE  WAVE  INVERSION 

Surface  wave  inversion  is  relativdy  straight  forward  once  the  dispersion 
is  correctly  defined.  This  requires  a  linearization  of  the  problem  and  the  sub¬ 
sequent  iterative  inversion  of  the  observations. 

Love-Wave  or  SH  Formulation. 

Fcdlowing  Keilis-Borok  et  al.  (1989),  the  equation  for  the  SH-wave  eigen¬ 
function,  in  a  cylindrical  coordinate  system  is 


d 

y(3) 

0  l/p. 

y(3) 

dz 

M 

pJk^-p«^  0 

where  k  is  the  wavenumber,  a>  is  the  angular  frequency,  p.  is  the  rigidity  and 
p  is  density.  The  boundary  conditions  for  surface  waves  are  that  p^z(z  =0) «  0, 
and  V^\z)-0,  as  2  -«>. 

A  discrete  solution  to  the  boundary  conditions  satisfies  the  equation 

Lnia^lQ-k%-l2  (12) 

0  0  0  V  1 

where  L  is  the  Lagrangian.  The  group  velocity  is  determined  from  the  rela¬ 
tion 


(13) 


If  the  medium  is  split  into  constant  velocity-density  layers,  of  thickness  d„, 
density  p„,  and  shear-wave  velocity  and  the  layer  is  bounded  by 
z  •  z„-d„  and  z  * z^/  Ihen  the  partials  of  the  phase  velocity  with  shear  velo¬ 
city  foe  fixed  fi’equency  and  density,  and  with  density  for  fixed  frequency  and 
shear-wave  velocity  are 


^a^Pm  *r 

W/o  ..zi-S, 


(14) 


and 


/  (V<^ydz  . 


(15) 
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JUyMfh  Vigoe  of  P-SV  Formulatiot. 

bi  a  cybndricai  coordinate  system,  the  equation  of  modem  governing  the 
fjgejiiirHactiotia  is 


0  l:X/(X+2|t)  l/(X+2»i)  0 

a  ^  0  0  1/M. 

dt  ftx  *  0  0  k  Pa 

p„  0  -pai*+4k*|4(X+|t)/(X+2|i)  -kX/(X+2p)  0 

whoc  the  vertical  component  eigenfunction  is  and  the  radial  component 
eigenfunction  is  V^.  The  dhptidty  at  the  free  surface  is  c  >  V'^\0)/V^^\0). 
The  boundary  cemditkms  few  surfaes  waves  are  that  p„{z >0)  >  0,  p^lz >0) »  0, 
V^\z)*0,  as  z  and  V^2)-0,  as  z  -«>. 

A  cUscrete  solution  to  the  boundary  conditions  satisfies  the  equation 


L-u%-k*li-2«2-./3, 
where  L  is  the  Lagrangian,  and 

(17) 

(18) 

Iff  [MV***)* + (X+2|iXV«)*]& 

(19) 

(20) 

(21) 

U 


The  group  vdk>dty  is  determined  frrom  the  rdation 


(22) 


If  the  medium  is  sfriit  into  constant  vdocity-density  layers,  of  thidmess  d„, 
(tonsity  p^,  compressional-wave  velodty  o^,  and  shear-wave  velocity 
and  the  layer  is  bounded  by  z  •*'<*  *  •  *w/  f*'*  partials  of  the 

phase  velodty  with  shear  vdodty  for  fixed  frequency,  compressions!  velodty 
and  density,  with  amipresdonal  vdodty  frxr  fixed  frrequency,  shear  velodty 
and  demity,  and  with  density  fcM'  fixed  frequency  and  compresdonal  and 
dicar  vdod^  are 


k  dz 


(23) 
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fi£.l  .  *f 

ldpj„  [  UIo  k  dz  }  k  dz 


[(V(Y  +  (V'^¥]<*2  •  (25) 


If  the  material  properties  within  each  layer  are  constant,  propagator  tech> 
niques  can  be  used  to  determine  the  dispersion  and  eigenfunctions  and  the 
integrals  can  be  anal3rtically  evaluated  (ftokrider,  1979;  Wang,  1981). 

Numerical  Partial  Derivation  Computaiton  and  Causality 

The  partial  of  group  velocity  with  respect  to  a  layer  parameter  v  is 


Co  9V  Co  Ic  df  dv  dT  dV  _ 


The  phase  velocity  partial  with  respect  to  period  is  numerically  computed 
using  the  relation 

3  dCo  (  9Co  9Co  1 

— - -  »  T  (— )t+|it-(— )t  at/21iT 

dT  dv  C  '  dv  ^  J 

In  these  expressions,  the  parameter  v  can  take  on  the  values  a,  ^  or  p,  and 
the  subscript  o  represents  the  value  in  the  parameter  in  the  purely  elastic 
model. 

If  a  causal  Q  is  introduced,  a  Futterman  causality  tied  to  a  reference  angu> 
lar  frequency  the  elastic  parameters  will  equal  the  anelastic  values  when 
the  angular  frequency  u  equals  .  The  causal  phase  velocity  c  is  given  by 

and  the  value  of  the  spatial  anelastic  attenuation  factor,  y,  is 

PB) 

In  order  to  obtain  the  correct  partial  derivatives  of  causal  phase  velocity  and 
gamma  with  respect  to  the  model  parameters,  the  following  are  used.  A  par¬ 
tial  or  a  parameter  without  the  subscript  o  represents  the  causal  value.  Note 
that  it  is  assumed  that  Qp/Qa  fe  «  fixed  ratio,  hence  a  partial  with  respect  to 
Qp*  may  involve  a  partial  of  with  respect  to  a. 


dc 

1  j 

dv 

ar  ^ 

dc 

.  JLwJjLi 

»  — m( — 
ir  tOf. 

[dp 

— ) 
“r 
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*1.  ,  -SL.  !^  q-i  .2-^— 
d!0»  26^  do  /.a. 


ta 

|dg 

c«  do 


ET  -  t4  |l  +  (2- 
oidhere  it  is  assumed  in  die  denvation  that 


Ce  (  Co  j  irw  j 

l^l’' 


is  negBgilte. 


dU 

dtZp 

€-Co 

do  ’ 

do 

Co 

[  Co 

r  ~  1 

dCp  Up 

do  c- 


t 


'irw 


Co  ‘"Qe 


_  1  .  ,  «  .  fc-Cp  V 

Co  I  J  ^ 


dU 


-1 


^®./5  ^9  %  dc 

‘^’cp  ^aOp* 


+  ii:^(if£.p  +  iStaOt) 

--‘'•p”  »a  0.' 


ir  Cm 


oi) 

(32) 

(33) 


(34) 


(35) 


An  iterative  linear  inversion  is  performed  because  of  the  non-linear  nature 
of  the  (MToblem.  At  any  stage  there  is  a  current  model  which  is  used  to  predict 
the  observations  and  there  is  also  a  lack  of  fit.  For  simplicity,  the  compres- 
^nal  wave  velocity  a  is  not  directly  inverted,  rather  it  is  tied  to  the  shear 
velocity  by  an  a  priori  assumption,  such  as  a  constraint  that  constant  Poisson's 
ratio  be  constant.  We  may  ^us  express  the  did^nce  between  observed  and 
predicted  values  as  a  linear  function  of  changes  in  shear  velocity  model  and 
changes  in  the  inverse  Q  model.  Since  the  e^ect  of  compressional  wave  Qp 
may  not  be  negUgible,  the  compressional  wave  Qp  is  related  to  the  shear 
^ve  (2p  by  a  ratio  Qp/Qp  for  the  layers.  Given  these  assumptions,  the 
inverskm  can  take  on  two  appearances,  non-causal  and  causal,  depending 
upon  whk^  partial  derivatives  are  used.  To  keep  the  notation  general,  the 
partial  derivatives  with  respect  to  layer  velocity  in  the  next  three  equations 
can  either  be  causal  or  non-causal. 

The  diiSdrences  in  observed  and  predicted  phase  velocities  are  modeled  as 


C^  ""Cffti 


Jc_ 

dg, 


dc 


*■  ■■  '  («) 


dc 


.-1 


^Qp. 
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the  difference  between  observed  and  predicted  group  velocities  is  given  by 


^obt  "“^pred 


dU 

9&t 


APi+ 


dU 


-1 


-  (37) 


3QS 


and  the  difference  between  observed  and  gamma  values  predicted  by  the 
current  model  is 


yobs’-lprtd 


-^3Lap,+  •  .  ■  +-^AP.  + 

ap,  ap,  aQg; 


•  (38) 


Once  set  up,  the  system  of  linear  equatioiw  is  solved  using  singular 
value  decomposition  (Russell,  1987).  Provisions  are  made  to  damp  the  solu¬ 
tions  as  a  mean  of  controlling  convergence.  Inversion  can  be  performed  in  an 
unweighted  manner,  or  can  be  perform  using  a  priori  weights  available  from 
the  error  estimates  from  the  dispersion  analysis.  In  addition,  a  differential 
inverse  can  be  use  to  solve  for  the  change  in  inversion  parameter  across  a 
layer  boundary,  rather  than  for  the  parameter  itself.  This  constraint  guaran¬ 
tees  a  smooth  model  and  avoids  unrealistic  alternating  low  velocity  layers  in 
the  resulting  model. 

MAINE  REFRACTION  EXPERIMENT 
In  1984,  the  U.  S.  Geological  Survey  ran  a  series  of  refraction  experi¬ 
ments  in  the  state  of  Maine  (Murphy  and  Luetgert,  1986;  Murphy  and  Luet- 
gert,  1987).  The  refraction  profiles  went  both  along  (NE)  and  cross  (NW)  the 
strike  of  the  Appaladtian  Mountains.  Surface  geology  varied  significantly  in 
lithology  and  geologic  age.  Because  of  the  excellent  data  set,  both  in  spatial 
sampling  and  calibrated  instruments,  this  data  set  was  analyzed  to  learn 
about  shallow  shear-wave  velocity  and  Q. 

Figure  1  shows  the  data  from  one  shot  point,  SPl  with  distance  increasing 
in  a  southeast  direction.  Fundamental  mode  phase  velocity  and  spatial 
attenuation  values  were  inverted,  with  the  resulting  earth  model  given  in 
Table  1.  To  avoid  any  bias  of  the  initial  model  upon  the  inverted  m^el,  the 
inversion  started  with  a  twelve  layered  model  with  velocities  and  densities  of 
each  layer  equal  to  that  of  the  halfspace.  The  inferred  vdodties  in  the  lower 
part  of  the  model  are  not  well  resdved.  The  high  values  are  not  unexpected 
given  the  highly  metamorphosed  geology  sampled  by  the  surface  waves. 

Figure  2  presents  the  fundamental  mode  synthetic  seismograms  for  an 
explosion  source  in  this  model.  The  arrival  times  and  waveforms  compare 
well  to  the  observed  data  of  Figure  1.  As  a  check  on  the  correctness  of  the 
inversion,  both  the  synthetic  and  observed  data  were  filtered  through  the 
same  narrow  band  band-pass  filters,  and  the  peak  amplitudes  plotted  as  a 
function  of  distance.  This  exercise  showed  that  the  derived  Qp  model 
explained  the  observed  decrease  of  amplitude  with  distance. 

FUTURE  DIRECTIONS 

Experience  with  limited  data  sets  indicate  certain  rules  for  estimating  the 
depth  limits  of  surface-wave  inversion.  Since  longer  periods  or  longer 
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DIST  (KM) 


Figure  1,  Date  ftw  SPl,  plotted  uwg  a  reduction  v^odty  of  6.0  km/a.  The 
date  have  been  low  pass  filtered  between  0.5  and  5.0  Hz. 


TABLE  1.  Vdodtv Modelfor SPl 


wnm 

iSSI 

mi 

fkni/s) 

(g/cmh 

Qa 

Qa 

0.25 

4.45 

2.57 

2.37 

64 

32 

0.25 

5.61 

3.24 

2.62 

64 

32 

0.25 

5.85 

3.38 

2.67 

64 

32 

0.25 

6.15 

3.55 

2.74 

70 

35 

0.25 

6.30 

3.64 

2.79 

74 

37 

0.25 

6.40 

3.70 

2.82 

84 

42 

0.25 

6.41 

3.70 

2.82 

96 

48 

0.25 

6.44 

3.72 

2.83 

114 

57 

0.50 

6.45 

3.72 

2.83 

142 

n 

0.50 

6.49 

3.75 

2.85 

190 

95 

0.50 

6.50 

3.76 

2.85 

284 

142 

6.65 

3.84 

2.89 

568 

284 

wavelength  attrface  waves  penetrate  deeper,  and  ^ce  loi^^  periods  are 
better  reaotved  at  laiger  distances,  two  such  rules  are  1)  the  surface  wave  pro¬ 
vides  leaolvable  infomation  to  depths  on  the  order  of  1/3  - 1  X,  where  X  is 
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Figure  2.  S3mithetic  traces  for  SPl,  plotted  using  a  reduction  vdodty  of  6.0 
kWs.  The  data  have  been  low  pass  filtered  betw^  0.5  and  5.0  Hz  and  also 
passed  through  a  2  Hz  seismometer  for  proper  comparison  with  Figure  1. 

the  longest  observed  wavdength,  and  2)  the  maximum  depth  is  on  the  order 
of  1/20  - 1/10  where  is  the  of  the  maximum  distance  of  observa¬ 
tion.  This  rdativdy  limited  depth  resolution,  compared  to  refraction 
analyds,  is  mitigated  by  the  ability  to  resolve  dumges  in  the  surface-wave 
waveguide  within  a  given  shot  spread. 

The  data  sets  andyzed  and  die  inverdon  programs  devdoped  both  con- 
ddered  land  based  data.  Introduction  of  a  water  layer  is  simple,  and  is 
already  encompassed  in  the  theory  presented.  There  is  no  difficulty  in  using 
surface-wave  inverdon  techniques  for  data  acquired  in  shallow  marine  sedi¬ 
ments. 
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ABSTRACT 

The  objective  of  this  study  is  to  demonstrate  that  a  laterally  varying  shal¬ 
low  shear-wave  structure,  derived  from  the  dispersion  of  the  ground  roll,  can 
explain  observed  lateral  variations  in  the  direct  shear-wave  arrival,  which 
occur  at  same  places  where  ground-roll  moveout  changes.  The  data  set  con¬ 
sists  of  multichannel  seismic  refraction  data  from  a  USGS-GSC  survey  in  the 
state  of  Maine  and  the  Province  of  Quebec.  These  data  exhibit  significant 
lateral  changes  in  the  moveout  of  the  ground-roll  as  well  as  the  S-wave  first 
arrivals.  A  sequence  of  surface-wave  processing  steps  are  used  to  obtain  a 
final  laterally  varying  shear-ivave  velodty  model.  These  steps  indude  visual 
examination  of  the  the  data,  staddng,  wavefcmn  inversion  of  selected  traces, 
phase  vdodty  adjustment  by  oross-cmrrdation,  and  phase  vdodty  inversion. 
These  models  are  used  to  ]»edict  the  S-nvave  first  arrival  by  using  2-D  ray 
tracing  techniques.  Observed  and  calculated  S-wave  arrival  matdi  well  over 
30  km  long  data  paths,  where  lateral  variations  in  the  shear-wave  velodty  in 
the  upper  1-2  km  is  as  nuidi  as  ±8%  .  The  modeled  corrdation  between  tihe 
lateral  variations  in  the  ground-toll  and  S-wave  arrival  demonstrates  that  the 
problem  of  shear-wave  statics  in  explication  seismology  could  be  posed  as 
problem  of  surface-wave  propagation. 

INTRODUCTION 

The  observed  quality  of  seismic  signals  can  be  degraded  by  their  propa¬ 
gation  through  shallow-inhomogeneities  and  lateral  variations  in  the  earth's 
near-surface  structure.  Near-surface  structure  introduces  delays  in  the 
recorded  travd  times  which  can  cause  dgnal  misalignment  that  <totorts  the 
apparent  earth  structure  and  stratigraphy  (Wiggins  et  al,  1976;  Rothman, 
1986).  Since  near-surface  time  anon^es  can  have  an  adverse  effect  on  all 
subsequent  processing  and  interpretation  of  seismic  data,  an  independent 
constraint  on  the  sh^ow  shear-wave  vdodty  could  improve  the  seismic 
exploration  results.  Usually  this  is  done  by  a  statics  analysis  applied  to  the  S- 
wave  first  arrivals.  Another  possibility  is  to  make  use  of  the  ground  roll,  or 
surface  wave,  to  define  a  riiallow  shear-wave  velocity  model. 

It  is  well  known  from  earthquake  seismology  that  the  inversion  of 
surface-wave  dispersion  can  be  usc^  to  define  a  ^ear-wave  velocity  struc¬ 
ture.  One  can  use  either  Love  or  Rayleigh  wave  fundamental  and  higher 
mode  dispersion.  Because  of  the  large  distances  between  seismograph  sta¬ 
tions  in  outhquake  seisnudo^,  inverted  earth  modds  represent  padt  ava- 
i^es  between  seismograph  stations  and  the  inverted  model  thickness  is 
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vmtMf  tigniBcagi^  tew  than  the  station  sepftkm.  Mui  (19B4)  used  Love 
wave  ^Safenkm  to  determine  static  cocrectiong.  He  used  all  the  traces  in  the 
shot  epreed  to  define  a  disperrion  oirve  and  correspOTidtog  vdodtjr  modd 
averafed  over  die  sfnead  lengdi.  For  purpose  of  static  correctkms,  die  vdo- 
efty  flMKM  lor  each  shot  spread  was  assig^ied  to  a  horizontal  position  at  the 
center  of  the  q^read.  Thus  a  apadaUy  smoothed  statics  model  is  derived.  In 
order  to  iaqprove  the  static  correcthm,  less  ^fiatial  averaging  is  re^iired. 
Efforts  have  been  made  in  diis  direction.  Nolrt  (1990)  combined  waveform 
modMng  tediidques  with  a  partitioned  inversion  schmne  to  trace  lateral  vari- 
atkms  in  earth  moddl  beneath  a  NARS  experiment  in  Europe.  His  technique 
is  vaHd  for  smiA  lateral  variations  in  vdodty  model. 

A  refraction  survey  in  the  state  Maine  conducted  by  the  U.  S.  Geologi¬ 
cal  Survey  and  the  Gcdogical  Survey  of  Canada  crossed  complex  geological 
structure  diet  ykdded  observable  lateral  variation  in  die  moveout  of  the  P- 
and  Swave  first  arrivals  and  also  die  surface  wave.  These  data  provide  an 
<^»portunlty  to  see  if  it  is  possible  to  use  data  frmn  a  sfri^  shot  q»ead  to 
define  a  laterally  varying,  shalfow  shear-wave  v^Kity  modd  from  die 
surfacewave  dispersion  and^then  to  see  if  this  velocity  model  can  ea^lain  the 
<d)served  lateral  variation  in  m  shear-wave  first  arrivds. 

DATA  AND  ANALYSIS 

DataSet 

Murphy  arul  Luetgert  (1966,  1967)  describe  the  experiment  and  data 
acquisition.  Scientific  results  are  ^ven  in  Klemperer  and  Luetgert  (1967)  and 
Lue^ert  et  al  (1969).  For  our  paper,  two  data  sets  were  chosen:  one  from  diot 
SP4  in  a  southeast  direetkm,  whidi  is  denoted  SP4-POS,  since  it  is  in  the 
direction  SPS,  and  the  other  frmn  SPl  in  a  northwest  direction,  vdikh  we 
denote  as  ^1-NEG.  Rgure  1  shows  die  shots  aikl  the  trend  (rf  the  recording 
poettkms  frmn  them. 

Data  PfoctttiHg 

Rgure  2  presents  the  observed  and  modded  data  for  the  SP4-POS 
qncad.  Tlw  vertied  emnponent  time  histories  are  available  at  approdmatdy 
1  km  ^wdngs  and  are  uWd  out  to  34  km  for  this  study.  The  traces  are  low 
pass  filtered  at  5.0  Hz  with  a  two-pole  Butterworth  cau^  filter  to  emphasize 
die  surfioe  wave.  A  reduetkm  vekteity  oi  3.5  km/s  is  used  to  focus  attended 
on  the  S-wave  first  arrival.  It  is  obvious  that  there  are  laterd  cherries  in  the 
ntoveout  of  the  S-wave  first  arrivd  and  doser  inspection  shows  that  this  is 
also  true  for  the  surfsce-wave  arrivd,  the  fundamentd  mode  Rayld^  wave. 
To  focus  <m  the  surfsce  wave  arrival,  these  same  data  are  plotted  in  Figure  3a 
widi  a  reduetton  vcfodty  of  2.5  km/s.  The  laterd  variations  in  the  surfsce- 
wave  moveout  are  apparent. 

To  define  die  surface-wave  dispersion  where  there  is  no  evidence  of 
laterd  vnriatkm  in  the  disperskm,  one  can  process  the  data  with  a  p-w  stadc 
^lIcMedian  and  YedUn,  1961).  This  will  lead  to  a  disperdon  curve  that  can 
be  inverted  kx  an  ea^  inodd.  When  there  are  sharp  dianges  in  die 
wav^itide,  the  pw  stack  may  yidd  two  dispersion  curves  (Mc^tar  et  d. 
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Rg.  1.  Map  of  the  ttudy  region,  showing  the  shot  points,  SPl  and  SP4,  and 
the  data  spreads  from  these  riiot  p<^t8. 

1968).  If  the  changes  are  gradual  the  p-«  stadc  dispersion  curves  will  be  less 
wdl  deOned  and  an  average  dispersion  curve  can  ^  inverted  Sot  a  first  order 
earth  model  which  can  dien  be  used  to  define  an  average  phase-vdodty 
dispersion  curve  tor  the  next  processing  step. 

If  the  signal  is  assumed  to  consist  <mly  of  a  fundamental-mode  surface 
wave,  then  ^  rignal  can  be  analytically  dispersed  to  another  distance  by  a 
simple  frequency  dooudn  operatCMr.  Let  the  observed  signal  at  distance  r  and 
frequency  /  be  5(r,/).  AIm  let  the  wavenumber  and  anelastic  attenuation 
coefficients  be  defin^  by  k{f)  and  y(f).  Then  the  signal  to  be  expected  at  a 
distance  f]  is 

If  the  predicted  time  history  at  distance  f]  is  plotted  at  r,  then  any  variaticm 
in  the  wavefcxrm  would  rquesent  a  vk^tkm  of  the  assumptions  of  wave  pro- 
{Uigatkm  in  a  laterally  unifcmn  medium  and  of  knowing  Ute  true  phase  vdo- 
dty  dispersion.  Figiure  3b  shows  the  result  ot  applying  such  an  operation  to 
the  trMBS  oS  Figure  3a  when  a  refereroe  distance  rj  •  30  km  is  used.  Since 
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Fig.  2.  OlHwrved  (light)  and  predicted  (heavy)  setsmograma  of  the  SP4>POS 
profile.  Both  trroea  have  bem  bandpass  filtered  betv^n  0.5  •  5.0  Hz.  The 
predicted  S-^ve  arrival  tunes  are  indicated  by  the  shaded  curve. 

die  Artt  Older  efiects  of  dtoperskm  have  been  equalized,  die  wavefonns  are 
slin&ar.  Differences  at  short  distances  are  due  to  lack  of  separation  of  die 
surface-wave  from  die  rest  the  signal  in  the  original  data.  At  large  dis¬ 
tances,  h%}i  frequency  noise  is  introduced  because  the  attenuation  operator 
must  reintroduce  high  frequencies  to  the  signal  in  equalizing  the  signal  to  the 
shorter  leferehce  distance.  Since  the  pulse  shapes  are  relativriy  independent 
(rf  distance,  the  first  order  effects  of  dtopersion  have  been  accounted  for.  Hg- 
lire  3b  can  be  vtewed  as  like  the  r^uced  travri  time  fdot  of  refraction 
setenology,  adiere  rdative  ofrsets  of  arrivals  represent  faster  or  sloiver 
arrivals.  It  is  obvious  from  Figure  3b,  and  also  from  Figure  3a,  that  diere  are 
dianges  in  the  surface-wave  moveout  at  about  12,  20  and  27  km  from  the 
source.  By  omqiaring  receiver  locations  to  gerdogical  maps,  these  dianges 
were  fourid  to  corrdate  to  mapped  dianges  in  surface  geology. 

The  initial  farocessing  steps  are  described  by  Herrmann  and  Al-Eqal» 
(1991).  Mult^de  filter  analysis  is  performed  on  eai^  trace  in  carder  to  estimate 
the  group  vdkidty  and  spectral  amplitude  of  the  fundamental  mode  as  a 
functfcm  of  frequiticy.  The  ^lectral  amplitude  data  are  then  used  to  estimate 
die  andarik  attenuatton  coefficient  y  that  appears  in  ecpiation  (1).  A  p-«» 
dack  is  perfornwd  to  obtain  phase  velodty  as  a  function  of  frequency.  These 
data  can  be  inverted  for  shear  vdodty,  ^  and  the  shear-wave  cpiality  factor. 
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Fig.  3.  (a)  Filtered  observed  data  of  Figure  1,  plotted  with  a  reduction 
velodty  of  2.5  km/s  to  focus  on  the  surface-wave,  (b)  Traces  of  Figure  3a 
equalized  for  dispersion  to  a  reference  distance  of  30  km.  This  display 
diminates  the  distraction  of  the  dispersion  in  order  to  focus  on  the  laterd 
variation  apparent  in  the  data  set. 

Qp  as  a  function  of  depth.  If  lateral  variations  in  dispersion  are  observed  in 
the  data,  the  spread  is  partitioned  into  segments  of  uniform  dispersion,  based 
on  an  examination  of  data  presentations  such  as  those  of  Hgure  3.  Waveform 
inversion  (Nolet,  1990)  is  applied  to  the  traces  at  the  segment  boundaries, 
and  the  resultant  modds  are  used  to  generate  fundamentd  mode  synthetics, 
which  are  then  cross  correlated  using  the  dispersion  curves  from  the 
waveform  matching  models  to  guide  the  intertrace  phase  unwrapping.  The 
result  is  a  dispersion  curve  spe^c  to  each  segment.  The  cross  corrdation  of 
synthetics  ra^er  than  actual  data  traces  3adds  smoother,  better  behaved 
intertrace  disperdon  curves.  These  dispersion  curves  are  independently 
inverted  to  provide  the  vertical  vdodty  distribution  for  each  segment. 
Finally  synthetic  seismograms,  induding  instrument  response,  are  generated 
for  the  resultant  laterally  varying  earth  modd  under  the  assumption  of  no 
mode  conversion  at  segment  bowdaries  (Keilis-Borok  et  d,  1989)  and  com¬ 
pared  to  the  observed  data.  A  laterally  varying  modd  is  judg^  appropriate  is 
the  predicted  and  observed  wavefemns  agree  both  in  shape  and  levd. 

Inoenknu 
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Table  1  riiows  the  remiltant  earth  model  for  the  SP4  profile.  In  this 
table,  eadi  segment  is  defined  by  the  distance  of  its  boundaries  from  the  shot 
point  A  sin^  Q  model  wu  assumed  for  the  entire  path  because  of  lade  of 
qpatid  resdutkm  inherent  in  estimating  Q  in  terms  of  a  ^>atial  decrease  in 
ampfitude.  fat  achSition,  is  %vas  assumed  that  the  P>wave  velodty,  a,  is  givoi 
by  a  «  1.732fi  arui  that  ■  2Qp.  Figure  2  oennpares  the  observed  (light 
fines)  and  fundamental  mode  surface-wave  traces  (heavy  lines)  for  this  data 
set.  Botit  sets  of  traces  have  been  bandpass  filtered  betv^n  0.5  and  5.0  Hz. 
The  agreemait  is  relativdy  good,  except  that  the  surface-wave  model  does 
not  pradid  the  later  reverbmtions  seen  in  many  observed  traces. 
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Another  way  of  presenting  the  results  is  shown  in  Figure  4.  Here  an 
OMSui  vdodty  in  each  layer  is  obtained  and  the  percentage  deviation  from  the 
mean  in  each  segment  is  displayed.  A  zone  of  high  velodty  is  seen  in  the 
third  segment  and  it  is  also  seen  that  segments  5-7  have  lower  surface  i^lod- 
tics  than  segments  1-4. 

A  question  arises  as  to  how  much  of  the  resultant  velodty  model  is  actu¬ 
ally  required  or  even  justified.  At  shc»t  distances,  little  dispersion  can  be 
resolved,  and  there  should  be  less  faith  in  the  resultant  model.  The  4  km 
thidc  modd  was  an  initial  assumption  in  the  inversion,  which  started  with  a 
hallqMce  model.  The  dispersion  data  were  permitted  to  smoothly  change  the 
modd.  However,  if  the  ^persion  data  could  not  provide  infnrmtion  about 
deepm*  structure,  that  structure  was  not  changed  much.  While  this  approach 
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Fig.  4.  Representation  of  lateral  varying  velocity  model  for  the  SP4-POS  data 
set.  The  shading  indicates  the  percent  variation  of  the  layer  velocity  about 
the  average  shear-wave  velodty  for  that  layer.  The  segment  number  is 
identified  by  the  name  at  the  top. 

might  not  lead  to  correct  structure  at  depth,  it  did  permit  the  data  to  define  a 
model  from  an  unbiased  starting  mod^  e.g.,  a  halfrpace,  and  prevented 
artificial,  inversion  induced  low  velocity  zones. 

To  provide  further  insight  both  on  depth  resolution  and  also  to  test  the 
possibility  of  modeling  the  S-wave  first  arrivals,  the  segmented  model  was 
used  to  predict  first  arrival  times  through  ray  tracing.  We  used  terven^s 
SEIS81  package  to  do  this.  For  simplicity  the  segment  model  was  assumed  to 
represent  a  single  2-D  laterally  varying  vdodty  model.  The  first  and  last  seg¬ 
ments  were  assumed  to  extend  b^nd  the  data  range  laterally.  Horizontal 
node  positions  were  placed  at  the  midpoints  of  each  segment,  and  vortical 
node  positions  wore  placed  at  the  midpoints  of  the  first  6  layers,  at  the  sur¬ 
face  and  at  a  depth  of  3.0  km.  The  velocity  in  each  layer  and  segment  was 
assigned  to  the  respective  nodes,  except  that  the  vdodty  at  ^e  surface 
equaled  that  of  the  top  la]rer  of  Table  1,  and  the  velocity  at  a  depth  of  3.0  km 
equaled  the  velocity  of  layer  8  of  Table  1.  This  latter  choice  was  related  to  the 
ad  hoc  experience  from  su^ce-wave  inversions  that  the  maximtun  depth  reso- 
lutk>n  is  on  the  order  of  one-third  the  longest  wavelength  resolved  or  one- 
tenth  to  one-twentteth  of  the  spread  length. 

Hgure  5  shows  the  rays  so  generated.  Their  arrival  times  are  plotted  in 
n^^ure  2  as  the  shaded  curve.  It  is  seen  that  the  predicted  S-wave  first  arrivals 
matdi  die  otxierved  arrivals  in  time.  Recall  diat  the  S-wave  arrival  times  were 
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Rg.  5.  Result  oi  tracing  rays  through  the  velocity  model  of  TaUe  1. 

not  used  to  derive  the  laterally  varying  vdodty  modri. 

Figure  6  shows  the  data  set  from  the  SPl  spread  to  the  northwest.  The 
shear-wave  arrivals  are  wdk  devdoped  and  show  lateral  variations  in  first 
arrival  times.  The  results  of  the  surface-wave  inversions  are  shown  in  Table  2. 
The  predicted  fundainentai-iiK)de  surfiice-wave  arrivals  are  also  shown  in  Fig¬ 
ure  6.  Rgure  7  shows  the  percentage  of  lateral  variation  in  each  layer  for  this 
nuxld.  Finally  Rgure  8  shows  the  rays  traced  through  the  model.  The 
predicted  S-wave  ^t-arrivals  are  shown  in  Figure  6  by  the  near-h<Hizontal 
shMied  curve. 

DISCUSSION 

This  study  focused  on  the  evaluation  of  the  shallow  shear-wave  vdodty 
structure  by  using  high  frequency  Rayleigh  surface  waves  recorded  as  part  of 
a  r^racthm  mqperiment.  Laterally  varying  earth  models  were  requited  to  fit 
the  data.  Seven  s^;ments  were  required  to  explain  the  latnal  variations 
observed  akmg  die  SP4-POS  (34  km  spread)  and  sbc  segments  required  for 
SPl-NEG  ^  km  8{»ead).  The  number  of  the  required  segments  is  a  function 
of  die  geologic  complexity  and  the  amount  of  change  in  vdodty  in  die  hor- 
iumial  directkm.  Most  oi  the  derived  lateral  changes  correlate  strongly  with 
puUttdied  gecdogic  data,  such  as  the  bedrodc  geological  map  of  Maine 
(CMberg  et  al,  1985).  GMierally,  observed  variations  in  velocity  widiin  the 
upper  1  km,  range  betuveen  ±1.2%  and  ±8.2%.  The  S-wave  first  arrivals 
show  ifignificant  lateral  changes  in  moveout  in  the  data  sets  at  places  where 


-20- 


11.00  10.00 
DIST  (KM) 


10.00 


14.00 


11.00 


11.00 


Fig.  6.  Observed  (light)  and  predicted  (heavy)  seismograms  of  the  SPl-NEG 
profile.  Both  traces  have  been  bandpass  filetered  between  0.5  -  5.0  Hz.  The 
predicted  S-wave  arrival  times  are  indicated  by  the  shaded  curve. 

strong  lateral  variations  have  been  observed  in  the  surface  wave  (Figures  2 
and  5).  Theoretical  S-wave  first  arrival  times,  as  obtained  by  2D  ray  tracing, 
show  strong  correlation  with  observed  S-wave  times  for  the  two  data  set  pro¬ 
cessed. 

In  the  case  of  the  SP4-POS  profile  we  note  that  Segment  1  has  the 
highest  velocity  and  consists  of  bands  of  Upper  and  Lower  Silurian  interbed- 
d^  pelite  and  limestone  and/or  dolostone.  Segment  2  consists  of  Devonian- 
Silurian  calcareous  sandstones.  Segments  3  and  4  have  alternating  bands  of 
Upper  and  Lower  Silurian  interbedded  pelite  and  limestone  and/or  dolostone. 
Se^ents  5  and  6  contain  alternating  bands  of  sulfidic  pelite,  Silurian  inter¬ 
bedded  pelite  and  limestone,  and  Silurian  interbedded  pelite  and  sandstone. 
Finally  Segment  7  has  Silurian-Ordovidan  calcareous  sandstone.  The  rela¬ 
tively  lower  velocities  in  the  upper  0.5  km  velocities  of  segments  6  and  7 
seem  to  correlate  with  the  sandstones.  Also  note  from  Figure  4  that  Seg¬ 
ments  1  and  3  have  similar  velocity  distributions  in  agreement  with  the  sur¬ 
face  geology. 

Along  the  SPl-NEG  path,  the  model  not  only  predicts  the  direct  shear- 
wave  arrival  and  Rg-wave  traces  that  match  the  observed  traces,  but  the  verti¬ 
cal  segment  contacts  agree  with  mapped  geologic  contacts  in  the  area  and 
also  with  the  interpretations  drawn  from  processing  reflection  data  recorded 
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Fig.  7.  Represmtation  of  lateral  varying  velocity  model  for  the  SPl-NEG  data 
set.  The  shading  indicates  the  percent  variation  of  the  layer  velocity  about 
the  avtai^e  shear^wave  velocity  for  that  layer.  The  segment  number  is 
identified  by  the  name  at  the  top. 

para&el  to  diis  path  (Luetgert  et  al,  1^9).  The  boundary  between  Segments  1 
and  2  maths  the  weston  boundary  of  the  Chain  Lake  I^ssif  and  the  position 
Western  boundary  Fault.  S^ment  2  omtains  Devonian  turbklites. 
The  Sandy  Stream  Fault  forms  a  boundary  between  Segments  2  and  3.  Seg¬ 
ment  3  cxHisists  of  the  Clinton  and  Chesham  formations  (volcanics,  sedi¬ 
ments,  mtiange).  Segments  4  and  5  contain  die  Frontenac  formatkm  (terri- 
geiKNis  and  volcanic  rodcs).  Finally  Segment  6  consists  of  the  highly  faulted 
and  fcMed  Silurian  and  Devonian  turbklites  of  Connecticut  Vaf^  -  Gasp4 
Syndhiorium  and  has  the  lowest  velocity  within  this  padi.  The  location  of 
ocmtacts  between  Segments  1  and  2  and  between  Segments  2  and  3  agree 
with  mapped  geological  contacts  and  with  events  C  and  Y  in  the  Luetgert  et 
al  (1909)  reflecfion  interpretation. 

The  ray  tracing  accomplished  two  tasks.  First,  it  demonstrated  the  power 
tA  die  surface-wave  based  shear-wave  velocity  model  to  predict  S-wave  first 
arrivals.  Second,  it  provided  some  idea  of  the  depth  of  penetration  of  the 
slmur-waves.  If  we  assume  that  the  first  arriving  shear  waves  penetrate  as 
deep  or  deeper  than  the  surface  waves.  This  conjecture  may  be  justified  by 
mode-ray  dtiality  considerations. 
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Fig.  8.  Result  of  tracing  rays  through  the  velocity  model  of  Table  2. 

While  the  data  analysis  satisfied  the  objective  of  this  paper,  additional 
work  is  required  to  make  this  a  robust  tool  for  exploration  geophysics,  espe¬ 
cially  with  respect  to  shear-wave  statics.  The  primary  problem  is  that  many 
individual  steps  are  required,  each  of  which  requires  significant  data  manipu¬ 
lation.  If  laterally  varying  phase-velocity  dispersion  functions  can  be  defined 
directly  from  the  trace  data  in  a  single  processing  step,  then  the  inversion  for 
the  corresponding  shear-wave  velocity  modd  is  relatively  trivial.  One 
approach  would  be  to  build  upon  the  combination  of  multi-channel  process¬ 
ing  and  phase  match  filtering  used  by  Russell  (1987)  for  anal3rzing  dispersion 
in  a  laterally  homogeneous  environment  to  model  the  interstation  phase 
dela3rs  using  a  laterally  varying  phase  delay  model.  Another  aspect  of  pro¬ 
cessing  required  to  malce  the  resiilts  more  convenient  would  be  to  keep  layer 
velocities  ^ed  and  to  invert  for  layer  thidcness.  This  would  lead  to  a  shal¬ 
low  shear-wave  velocity  model  better  related  to  statics  modeling.  Keilis- 
Borok  et  al  (1989)  provided  the  necessary  surface-wave  perturbation  theory  to 
accomplish  this.  Thus  the  routine  use  of  the  ground  roll  to  constrain  shallow 
shear-wave  velocities  for  use  in  shear-wave  statics  estimation  is  possible. 
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ABSTRACT 

Knovde^e  df  aeismic  shear  wave  vdodty  and  Q  is  required  fmr  both  seismic 
hazard  analysis  uui  the  dting  of  critical  structures  since  Utese  parameters 
gcwmm  die  teaiiS8daai<Mi  <ii  the  seismic  signal  from  the  earthquake  to  die  site 
and  also  control  the  site  response  itsdf .  These  parameters  vary  spatially  and 
may  not  be  readity  availaUe.  The  purpose  of  diis  paper  is  to  show  how  these 
can  be  estimated  using  surfrioe-wave  analysis  tediniques. 

Data  £rom  USGS  refraction  surveys  in  Maine  were  selected  which  show 
dispersed  short-period  Rg  waves  <fundam«(ital  mode  Rayleigh  waves).  Pro- 
cessiiig  tediniques  appli^  to  the  data  indude  filtering,  wai^cnm  inversion 
of  selected  trace,  phase  vdodty  stacking,  interactive  amplitude  processing, 
invecekm  of  disperrion  parameters  fmr  a  layered  shear-wave  vdodty  and  Q 
modd,  and  finafly  verifying  the  derived  vdodty-attenuation  modds  by  com¬ 
paring  bodi  shape,  ab^ute-amplitude  and  arrival  times  of  syndietic  and 
obsen^  time  series.  The  resultant  vdodty  modds  show  that  the  shear- 
wave  vdodty  varies  between  2.43-2.81  km/sec  for  surface  layer  and  Qp  is  less 
than  50  in  the  upper  kilometer  of  the  metamorphic  terrain  of  Maine.  The 
modds  oorrdate  with  the  type  of  geologic  rocks  encountered.  While  the 
study  focused  on  upper  crustal  structure,  the  techniques  and  eqierimental 
prooi^ure  can  be  used  for  shallow  rite  characterization  by  appropriate  scaling 
of  the  data  aoquidtion. 


INTRODUCTION 

The  analysis  of  surface  wave  dispersion  provides  a  convenient  method 
fm  determining  the  elastic  properties  of  a  me^um.  One  can  use  dther  Ray¬ 
leigh  or  Love  wave  signals.  By  inverting  their  phase  and  group  vdodty 
dispmion  and  firequency  dependent  anelastic  attenuation,  a  shear-wave  velo- 
dty  and  Q  model  can  be  obtained.  The  depth  extent  and  vertical  resolution 
of  the  earth  models  is  a  function  of  the  largest  and  smallest  frequendes 
observed  in  the  signal.  Shallow  refraction  surveys  and  quarry  tdasts  are 
good  generators  of  Rg  waves,  the  short  period  Rayldgh  wave.  These  waves 
are  characterized  by  their  huge  amplitudes  and  dear  dispersion.  The  Rg 
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waves,  because  of  th^  high  frequencies  relative  to  long  period  teleseismic 
observations,  are  confined  to  sh^ow  depths,  which  makes  these  waves  a 
very  good  diagnostic  tool  in  studying  the  shdlow  structure.  Typically  they 
will  not  sample  earth  structure  as  deeply  as  the  observed  P-  and  ^wave 
arrivals. 


THE  SEISMIC  REFRACTION  DATA  SET 

In  the  fall  of  1984,  the  USGS  collected  850  km  of  seismic  refraction  data 
in  Maine  (Murphy  and  Luetgert,  1986;  Murphy  and  Luetgert,  1987).  It 
includes  one  major  profile  of  more  than  300  km,  running  perpendicular  to  the 
strike  of  the  dominant  structure  and  crossing  the  major  tectono-stratigraphic 
units  of  east  and  central  Maine,  and  four  shorter  profiles  extending  parallel  to 
the  strike  and  sampling  individual  tectonic  units.  Each  profile  consists  of 
several  observation  lines.  Each  line  induded  several  shots  with  the  120  verti¬ 
cal  component  instruments  deployed  spatially  over  90  km  spreads.  The 
receiver  and  shot  locations  are  t^eved  to  be  accurate  within  15  m.  The  error 
in  geophone  location,  relative  elevation  and  timing  are  considered  negligible 
and  less  than  the  error  in  travel-time  due  to  the  changes  in  weathering  zone 
velodty  or  thickness  (Klemperer  and  Luetgert,  1987).  Two  data  sets  are 
selected  to  illustrate  the  use  of  surface  waves  in  determining  shallow  shear- 
wave  velodty  and  Q  structure.  These  are  from  profiles  SP13-NEG  and  SP26- 
NEG,  which  are  parallel  to  the  structural  trend  in  Maine.  Figure  1  shows  the 
two  shots  and  their  observations  p<^t8.  These  two  profiles  were  chosen 
because  of  the  lack  of  observable  lateral  variation  in  the  surface  dispendon 
along  the  paths.  These  simple  data  sets  are  displayed  in  Figures  2  and  3.  The 
obs«^ed  data  have  been  bandpass  filtered  between  0.5  and  5.0  Hz  using  a 
four  pole,  two  zero,  Butterwmrth  filter.  The  data  are  plotted  as  a  function  of 
reduced  travd  time  and  distance.  For  the  SP13-NEG  profile,  distinct  P-  and 
S-wave  arrivals  can  be  seen. 

REGIONAL  GEOLOGIC  SETTING 

Most  of  the  State  of  Maine  is  underlain  by  tighdy  folded  metasedimen¬ 
tary  sequences  of  Ordovician  to  Devonian  age.  The  major  structural  feature  in 
which  the  SP13-NEG  path  falls  is  the  Central  Merrimack  Syncdinorium.  The 
syndincnrium  consists  of  tightly  folded  faulted  metasedimentary  secjuences 
and  turbidites  of  Ordovician  to  Devonian  age.  The  SP26-NEG  path  falls 
within  the  volcanic  bdt  along  the  coast  of  Maine.  The  surface  geology  along 
the  profiles  varies  significantly  in  lithology  and  geologic  age  reflecting  the 
complicated  geologic  past.  The  rocks  of  SP13-NEG  path  are  interbedded  pd- 
ite  and  sandstone,  whereas  those  of  the  SP26-NEG  are  mostly  voka^cs 
(Osberg  and  others,  1985). 


-71.50  -69.50  -67.50 


Hf.  1.  thowrlng  Om  Stale  of  Maine  and  the  locatkma  of  the  two  shot 
poiiila  Ortangles)  and  tiie  observation  points  for  eadt  shot  point  Shot  point 
SPlS'^fEG  is  in  dte  souOiwcst  part  of  die  state,  while  SP26-NEG  is  in  die 
easleni  part  of  the  Stale. 

DATA  PRCXXSSING  AND  ANALYSIS 

The  raowded  multichannel  refraction  data  sets  show  excellent  Rg-wave 
diqpetakm.  Processing  tediniques  were  devdoped  diat  could  be  apjdkxl  effi- 
denlfy  to  obtain  die  desired  shear-wave  vdodty  and  Q  structure.  The  tech- 
nfcptos  are  based  on  either  peifonning  direct  inversion  of  waveforms  for  the 
sheer-wave  vdodty  modd  or  extracting  the  dispo^ion  and  anelasdc  attenua¬ 
tion  ooeflklents  that  are  then  inverted  for  the  shear-velocity  and  Q  modds. 
The  data  sets  are  first  p-w  stadeed  to  obtain  phase  velocity  curves.  The  ane- 
lastic  oodfident  data  are  obtained  after  extracting  spectrd  amplitudes  by 
applying  muldple-filter  andysis  followed  by  interaedvdy  processing  the 
aa^tudes  to  determine  their  spatia]  anelastic  attenuation.  Tlie  phase  vdo- 
dto  and  andtotk  dispersion  data  are  then  inverted  to  obtain  the  shear-wave 
vdodty  md  Q  mod^.  While  the  multiple  filter  analysis  also  yidds  useful 
group  vdodty  information,  this  is  not  used  because  of  the  extent  phase 
vdodty  data  and  also  because  of  the  desire  to  demonstrate  that  the  earth 
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Rg.  2.  Data  fat  uaad  lor  ^344EG  profila.  IKa  tracat  are  tfie 
obaarvad  data.  The  heavy  traces  are  the  forward  syntihetic  data  derived  from 
the  inverted  velocity  and  Q  modds. 

model  can  {nredkt  the  observed  time  histories. 
p  ~ta  Stacking 

The  p  staddng  technique  is  used  to  obtain  phase  vdodty  dispersion 
without  making  asmimptions  about  source  ph^  or  the  number  of 
wavdengths  between  the  source  and  single  recdver,  both  of  vddch  add  ambi¬ 
guity  to  one-  or  two-station  phase  velocity  estimation  techniques.  This  tech¬ 
nique  provides  excellent  results  vdien  processing  multidtannd  data  by  being 
aUe  to  isolate  individud  normd  modes  of  propagation.  It  invdves  two 
transformations,  a  p  -  r  stack  followed  by  a  Fourier  transform  over  r. 
(McMechan  and  YedUn,  1961;  Russell,  1967;  Mokhtar  and  others,  1968; 
Hemnaim  and  Aleqabi,  1991). 

The  seismic  waves  c^Merved  at  a  distance,  x„,  can  be  represented  as  a 
sum  of  superimposed  imxles  as 

m 
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Fig.  3.  Dili  Ml  uMd  for  tti*  SP26-NEG  profile.  The  U^t  trace*  axe  fiie 
oimreeJ  data.  The  heavy  trace*  ate  die  forward  syndielk  data  derived  from 
dte  inverted  velodty  and  Q  modeb. 

ttdicre 

A„(x,u)  ■  /(»)  S,(»)K.(»)  .-“-(-Xw  ,  p) 

m  la  tfra  mode  number,  /(w)  the  instrument  response  (complex),  S^(tA)  the 
source  spectrum  (curniplex),  the  path  response  (real),  and  7m(<A)  the 

atteitustion  ooeIRdent  (real).  The  wavenumber  k  is  rdated  to  the  phase  velo¬ 
city  c  by  the  reiatkm  k  >  m/c  .  If  the  Fourier  transform  of  (1)  is  defined  as 
the  p  -«»  stadc  is  defined  here  as 

Fo>.“)  -  2  (?) 

* 

or  ex|dicitiy  as 

•  22  ■*-*•**■.  (4) 

R  m 

If  the  s%nal  consists  of  a  single  mode,  then  it  is  ea^  to  see  that  this  summa¬ 
tion  becomes  extremal  if  up  ■  1^.  By  evaluating  I  F(p,<a)  I  as  a  function  of 
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p  kx  »  fixed  frequency,  one  can  search  for  the  values  of  p  that  make  the 
expression  extrenud.  p  is  related  to  the  phase  velocity,  c,  by  the  simple  rela¬ 
tion  p  >  1/c .  For  m^timode  signals,  secondary  extrema  will  occur,  which 
may  give  the  dispersion  of  another  mode.  Since  modal  ecdtation  is  frequency 
dependent,  different  modes  will  dominate  at  different  frequendM. 

Multiple  Filter  Analysis 

Multiple  filter  analysis  was  introduced  by  Dziewonski  and  odiers  (1969) 
for  the  purpose  of  extracting  the  surface-wave  group  velodties  and  spectral 
amplitudes  from  dispersed  waveform.  Define  the  Fourier  transform  of  a 
dispersed  propagating  seismic  signal  at  distance  x  by 
and,  a  Gaussian  fflter  as 


H(«) 


0 


I  I 
I  <D-&>o  I 


(5) 


where  ,  a  is  a 

filter  parametn  controlling  the  tradeoff  between  time  domain  and  frequency 
domain  resolution,  mq  is  the  filter  center  frequency  and  is  the  filter  cutoff 
frequency.  The  result  of  appl3ring  this  filter  to  the  dispersed  signal  is  a  time 
domain  signal  whose  envelope  reaches  a  maximum  at  the  group  velodty 
arrival  times  and  given  is  by  Herrmann  (1973)  as 


2ir  a 


(6) 


where  Uq  is  the  group  velodty. 


Computationally,  the  time  series  is  filtered  using  (5),  a  Hilbert  transform  of 
the  filtered  time  series  is  taken,  aiul  the  two  are  combined  to  form  the 
envelope.  The  envdope  is  searched  for  the  relative  maxima,  gouw'  which 
the  group  velodty  is  computed.  The  spectral  amplitude  is  estimated  from  (6) 
by  the  relation 

(x  ,b)o)  >  gin«x(2‘>r/«oKo/'ir)’‘ 


A  given  set  of  seismic  traces  can  be  processed  using  multiple  filter  analysis 
individually  or  as  a  group.  Assuming  a  laterally  homogeneous  waveguide, 
the  result  of  a  single  trace  processing  could  be  utilized  to  obtain  adequate 
earth  structure.  For  a  group  of  traces  we  could  use  two  approaches  for 
estimating  the  average  vdodty  and  anelastic  attenuation.  We  can  run  multi¬ 
ple  filter  analysis  fc»  each  trace  and  consolidate  the  individual  group  velodty 
and  spectral  amplitude  estimates  into  a  large  file  and  use  an  interactive  com¬ 
puter  program  to  manually  isdate  modes  on  the  basis  of  a  group  vdodty 
window,  and  then  use  all  information  associated  with  that  wi^ow  to  deter¬ 
mine  a  mean  group  vdodty  and  an  anelastic  attenuation  coeffident,  y,  from 
the  decay  of  amplitiKle  with  distance.  The  modd  use  to  obtain  y  is  that  the 
qpectrd  amplitude  ot  a  dngle  mode  and  frequency  vary  with  distance 
acccardii^  to  the  rdation 
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A(x,«)  •  Aoe^*/Vz 

The  Moond  wmy  of  analyzing  group  of  traces  for  average  group  velocity 
mid  avonge  andastic  attenuation  co^dents  is  to  im^ement  dte  group  velo¬ 
city  deck  oonoqpt  oi  Barker  (1968).  The  idea  to  to  obtain  narrow  band  pass 
Gaussian  filler  of  eadi  trace,  in  a  manner  of  (5),  and  then  to  stack  the 
envdopcs  oi  diflcrent  traces  for  a  given  frequency  after  transfcnnning  the 
envdo^  from  function  oi  time  to  function  of  group  velocity.  A  search  then 
is  performed  to  find  the  stack  maximum  amplitude  which  defines  the  group 
velociQr.  Ush^  tlds  informatkm,  the  spectral  am^ditudes  associated  wito  the 
correspondiog  group  arrival  time  can  be  extracted  from  the  cxriginal  filtered 
traces,  and  a  regresdon  perfumed  to  cibtain  y.  Spectral  amplitudes  can  be 
used  to  cstfoiate  the  andmtic  attenuation  coefficient. 

SurfKt-Vfmm  Dtspenion  Inversion 

An  inqpBdt  theoretical  rdationship  exists  between  modes,  phase  and 
group  vdodties,  and  the  structural  ea^  paranwters  sudi  as  shear  vdodty 
«id  density.  Ity  combining  linear  inverdon  theory  with  variatiofud  tech¬ 
niques  to  evaluate  necessary  partial  derivatives  with  respect  to  modd  parame¬ 
ters,  die  data  can  be  invoted  to  obtain  an  earth  modd  (  Russdl,  1967).  A 
lineiur  iterative  inversion  sdimne  can  be  used  to  minimize  the  difference 
between  cibaerved  and  predicted  phase  vdodty, 

^ 

group  vdodty, 

+  . +  . . .  + 
and  andastk  attenuation  coeffidents, 

^ 

The  resulting  system  of  linear  equations  can  be  solved  using  inversion  theory 
to  find  shear-wave  velodty  and  Q  structure.  The  system  of  the  entire  set  of 
observed  data  will  have  the  matrix  furm 

Y-AX  +  e  (10) 

whuce,  Y  fr  an  m  z  1  is  the  vectcnr  of  differences  between  thedreticd  disper- 
skm  dbita  and  observed  dispersion  data.  X  is  an  n  z  1  is  vector  of  unknown 
parameters,  cisan  mzlisa  vector  of  residuals,  and  the  kemd  matrix  A  is 
Mt  SI  z  n  matrix  collating  of  the  partid  derivatives  of  dispersion  parameters 
with  respect  to  starting  modd  parameters.  By  the  least  squares  prindple,  a 
set  oi  nmM  parameters,  X,  is  conddered  a  valid  sdution  if  it  wfil  minimize 
die  sum  of  the  squares  of  the  residuals,  c. 

To  simplify  the  minimization  process,  Lawson  and  Hanson  (1974)  fritro- 
duced  a  tjngMhw*  value  decomposition,  wldch  is  an  orthogond  trandbnnation 
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of  the  matrix  A  into  three  matrices  and  can  be  written  as  A  «  UAV^  where^  U 
is  an  m  X  m  orthogonal  matrix,  V  an  n  x  n  orthogonal  matrix,  and  A  an  m  x  n 
diagonal  matrix.  Siiuie  an  "orthogonal  matrix"  is  a  matrix  whose  transTOse  is 
also  its  inverse,  the  solution  vector  X  is  given  by  X«VA~^lrY  or 
equivalently  X  ■  (A^  A)“'A^Y. 

In  situations  where  the  inversion  problem  is  poorly  constrained,  at  least 
one  of  the  singular  values  \j  will  be  small  or  zero,  causing  an  eccessive  mag¬ 
nification  of  the  corresponding  vector  v^  .  This  instability  in  the  solution  can 
be  controlled  by  adding  the  constraint  of  keeping  the  size  of  a^X^X  small. 
The  corresponding  solution  vector  is 

X  *  (A^A  +  a^ir^A^Y  (11) 

where  a  is  the  damping  value  and  1  is  the  identity  matrix.  If  the  damping 
parameter  a  is  increased,  more  weight  is  imposed  on  minimizing  the  solution 
norm,  and  less  on  reducing  the  least-squares  residuals,  and  if  it  set  =  0  the 
problem  of  the  stochastic  least-squares  will  be  reduced  to  standard  least- 
squares  problem.  In  our  formulation  the  solution  vector  X  represents 
changes  in  the  shear  velocities  and  inverse  Q,  and  is  used  to  update  the  start¬ 
ing  model.  The  process  starts  again  with  this  new  model  and  is  repeated 
until  the  desired  level  of  convergence  is  obtained. 

Y/aveform  Inversion 

The  waveform  inversion  procedure  uses  waveform  data  to  find  the 
model  parameter,  which  minimizes  the  difference  between  the  observed 
time  histories,  s(Xj,t)  and  the  theoretical  seismograms  (  represents 

model  parameters  and  could  include  velocity,  density,  Q  structure,  source 
location  and  mechanism.  Generally,  a  priori  knowledge  and  some  parameter 
assumptions  are  required.  Mathematically,  the  minimization  problem  can  be 
written  as 

F(e)  .  I  uit.x.Q-S{t,x)  I  Ut  (12) 

¥  0 


wl^ere  if  »  the  norm,  T»  time  span,  and  F({)»  the  objective  function.  The 
solution  to  the  above  optimization  problem  is  to  find  ^e  vector  (  (i.e.,  the 
model)  that  minimizes  the  function  F(().  A  detailed  discussion  of  this 
scheme  is  presented  by  Nolet  and  others  (1986)  and  Nolet  (1987,  1990). 


Synthetic  Seismograms 

In  a  homogeneous  earth  model,  the  Fourier  spectrum  for  single  mode 
recorded  at  receiver  distance  r  and  depth  z  firom  the  a  source  at  depth  h  is 


-4m 

ti(r,<^,z  ,«>)-«  * 


e-*-  €V(a>,z)  W(«o,jil 

V2irlT  {TcUIoi*  {TcUIf^li^ 


(13) 


wha«  (s  is  the  angular  frequency,  k  is  the  wave  number,  c  is  the  phase  velo¬ 
city,  U  is  the  group  velocity,  and  /q  is  the  energy  integral,  c  is  either  a  real  or 
imaginary  number,  depending  upon  the  component  of  motion.  V(tt>,z)  is  the 
eigenfunction  sampled  at  the  receiver,  and  IV(a»,2)  is  a  source  mechanism 
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UlMMF  comtoiatipii  of  dg|en£^l»^o|ls  at  the  source.  The  term 
A  » ¥  peM  the  ainplihide  response.  The  four  terma  in  the  eque- 
req>ectivel)r  the  phase  term  from  the  Hanhel  function,  the 
1^149$  term  froq^  Har^  function,  t)^  sampting  of  the  wavefield 
#  tliiasfQehwi|»  eoA  eapipling.  of  the  excitation  at  the  source. 

WSMflQ  mOCFTY  MODiLS 

l^jipe  4  the  multichannel  phase  velocity  sta<h  for  SP13-NEG  path. 
'Pm.  sorted  with  e  twelve  layered  earth  modd  with  velocities  and 

ofi  eech  layer  equal  to  that  of  the  halfspace.  Table  1  presents  the 
aheer-weve  velodty  values  for  this  path.  The  velocities  in  the  lower 
fail  of  t|ie  model  are  not  well  resolved.  Alro  shown  in  the  same,  table  are 
Ih#  valime  for  this  path.  These  values  were  obtained  by  inverting  the  ane- 
lastic.  aftenualion  coefficients,^,  obtained  from  the  interactive  processing  of 
Urn  mectrel  amplitudes.  Table  2  gives  the  shear^waye  velocity  model  for  the 
ol^mfrmd  from  iq>ptying  the  waveform  inversion  technique  to 
taace  at  a  distance  of  15.13  from  the  source.  The  model  for  path 
(laW#  WW  pbtained  by  inverting  observed  anelastic  coeffident 
obtafried  fr^  multichannel _^up  vdodty  and  andastic  coefficient  stadking 
tedmiqaes  (Barker,  1968).  'Ae  shear-wave  vdodties  of  the  SP26-NEG  path 
art  l^mr  than  those  of  the  SP13-NEG  path  aitd  this  is  expected  since  the 
SP^b'^fEG  path  lies  within  the  metavolcanic  rodcs  of  the  coastal  lithotectonic 
unit,  whereas  the  SP13*NEG  path  lies  within  the  Central  Merrimadc  Syndino- 
rium  with  aaostty  interbedded  ptKte  and  sandstone  rocks. 

SYNTHETIC  SEISMOGRAMS 

f  mr  tha  laterally  homogeneous  models  derived  in  this  study  (Tables  1 
and  Hi  aamlhttic  ^  wava  saismioyams  wart  calculated  asaumin^  a  step 
souma  function.  Tha  syntha^  seismograms  were  pam^  tiuough  a 
aacamd  aidai  Buttmwaftit  with  a  low-cut  comer  frequency  at  2  Hz  to 
ap|W<mtinala  tim  aanaor  response  of  the  original  data,  and  tfren  both  observed 
and  syrti^tic  seismograms  were  band-pass  filtered  between  0.5-5.0  Hz  to 
tmphasiw  ihe  Rg  wave.  Figures  2  and  3  compare  the  synthetic  fundamental 
taada  aaiamograma  and  tHa  original  data  for  paths  spi3-N£G  and  sp^h|EG, 
ftapatti^^-  An  aactUent  fit  is  obtained  in  both  surface-waya  s^pa  and 
tima  at  a  given  dfrianca-  Other  diapiaya  imve  shown  an  excefi^t  fit  to 
Mtk  ampfitudes  as  a  function  of  diafance,  thus  supportiim  Qp  model 
loundr 

PlfiQJfiSION  AND  CONCLUSIONS 

This  study  locuaad  on  presenting  geophysical  techniquea  for  the  eyahm- 
tinn  of  the  shaBow  shear  velocity  structures,  utilizing  high  frequen<ty  fiaX* 
la^  wavaa  (fig)  reomd«d  by  tba  USGS  during  the  Maine  refraction  expe^ 
inenL  A  variety  of  gepphyaical  proemsing  techniques  were  appliad  tp  the 
data  aa^  among  tfrem,  filing,  phase  vel^ty  sucking,  wavefom  invoraipn 
ad  aelppled  traees,  inversion  of  phmc  vdodty  and  synthetic  seismogram  calcu- 
latRm. 
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Fig.  4.  Result  of  the  p-<a  stadc  for  the  traces  in  Figure  2.  The  stack 
vidues  are  contoured  as  a  function  of  phase  velocity  and  frequency. 
Large  amplitudes  are  seen  along  the  locus  of  the  actual  dispersion 
curve,  for  what  may  be  higher  modes,  and  for  stack  noise  due  to 
spatial  and  temporal  aliasing.  The  fundamental  mode  phase  velocity 
dispersion  curve  from  the  inverted  model  is  also  shown. 

The  stacking  techniques  require  some  manual  intervention  to  correctly 
define  the  dispersion  data  for  a  given  mode.  This  is  required  since  different 
modes  are  dominant  at  different  frequencies.  This  complicates  the  dispersion 
data  selection.  Other  complicating  factors  arise  from  the  medium  geometry, 
lateral  velocity  variations  and  spatial  attenuation.  When  no  lateral  variation  in 
the  earth  model  is  expected,  then  a  single  well  dispersed  trace  is  sufficient  to 
define  the  earth  model,  as  was  done  with  the  data  of  Figure  3.  To  estimate 
the  anelastic  attenuation  of  this  data  set,  the  group  velocity  stacking  concept 
of  Barker  (1988)  was  used  to  automatically  estimate  the  anelastic  attenuation 
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Table! 

Velocity  and  Qp  Model  for  SP13NEG 


Layer 

Thkdcness 

km 

Compressional 

velocity 

knVsec 

Shear 

velocity 

km/sec 

Density 

gm/cc 

Q« 

Qp 

0.250 

4.218 

2.435 

2.350 

63.0 

31.50 

0.250 

4.990 

2.881 

2.498 

65.0 

32.55 

0.250 

5.094 

2.941 

2.519 

73.0 

36.50 

0.250 

5.493 

3.233 

2.598 

82.0 

41.00 

0.250 

5.636 

3.254 

2.627 

92.0 

46.00 

0.250 

5.725 

3.305 

2.645 

106.0 

53.00 

0.250 

5.756 

3.323 

2.651 

124.0 

62.00 

0.250 

5.767 

3.330 

2.653 

148.0 

74.00 

0.500 

5.770 

3.331 

2.654 

185.0 

92.50 

0.500 

5.770 

3.331 

2.654 

247.0 

123.50 

0.500 

5.770 

3.331 

2.654 

371.0 

185.50 

5.770 

3.331 

2.654 

743.0 

371.50 

coefficients,  y.  The  success  in  matching  the  observed  time  histories  indicates 
that  it  may  be  possible  to  automatically  process  such  data  without  too  much 
analyst  intervention. 

In  spite  of  limitations  of  the  processing  techniques  used,  frequency  con¬ 
tent  of  ^e  data,  and  the  geologic  complexity  of  the  area,  die  analysis  pro¬ 
vided  important  information  that  can  characterize  the  subsurface  along  the 
processed  paths  in  Maine.  Higher  velocities  are  observed  (2.816  km/sec  in 
top  layer)  for  paths  where  metavolcanic  rodcs  are  observed  (SP26-NEG). 
Lower  velocities  were  observed  (2.435  km/sec  in  the  top  layer)  on  the  SP13- 
NEG  path,  where  the  rocks  are  mostly  interbedded  sandstone  and  pelite. 
The  values  of  are  less  than  50  in  the  top  1  km  for  both  paths.  Su(^  low 
values  also  have  also  been  found  by  Mokhtar  and  Herrmann  (1988)  for  the 
Arabian  shield.  The  exact  cause  for  such  low  Qp  values  cannot  be  explicitly 
attributed  to  changes  in  lithology  but  may  be  related  more  to  the  formation 
pore  geometry,  fluid  content,  and  cradc  density  all  acting  together. 

The  teduiiques  presented  here  can  be  applied  to  problems  of  other 
scales.  For  example  if  the  refraction  experiment  is  performed  over  a  3  km 
spread,  and  frequencies  of  10  -  30  Hz  are  used,  one  should  be  able  to  invert 
for  the  shear  %vave  velocity  and  Q  in  the  upper  300m.  Shorter  spreads  would 
focus  on  proportionately  shallower  depths.  The  use  of  surface  waves  does 
not  require  dHlling  and  logging,  is  essentially  non-invasive  and  results  can  be 
obtained  fairly  rapidly,  with  all  the  processing  performed  on  an  ordinary  PC. 
The  shear-wave  velocity  is  directly  related  to  the  shear  modulus  of  the 
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Table  2 

Velocity  and  Qp  Model  for  SP26-NEG 


Layer 

Thickness 

km 

Compressional 

velddty 

km/sec 

Shear 

velocity 

km/sec 

Density 

gm/cc 

Qa 

0.250 

4.877 

2.816 

2.507 

56.0 

28.00 

0.250 

5.845 

3.372 

2.681 

59.0 

29.50 

0.250 

5.912 

3.410 

2.681 

65.0 

32.50 

0.250 

5.997 

3.460 

2.735 

77.0 

38.50 

0.250 

6.163 

3.555 

2.785 

90.0 

45.00 

0.250 

6.387 

3.684 

2.852 

105.0 

52.50 

0.250 

6.491 

3.745 

2.859 

123.0 

61.50 

0.250 

6.696 

3.863 

2.905 

147.0 

73.50 

0.500 

6.845 

3.949 

2.939 

184.0 

92.00 

0.500 

6.859 

3.957 

2.943 

246.0 

123.00 

0.500 

6.867 

3.962 

2.945 

369.0 

184.50 

6.885 

3.972 

2.950 

738.0 

369.00 

material  and  is  thus  an  important  parameter  for  engineering  design. 
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ABSTRACT 

Sturface  wave  data  from  eight  chemical  explosions  along  a  USGS  refrac¬ 
tion  profile  in  the  State  of  Maine  are  used  to  estimate  isotropic  moment  - 
yield  relations.  The  moment  estimates  are  obtained  by  using  the  surface-wave 
dispersion  data  to  define  an  earth  model,  the  decrease  of  surface-wave  ampli¬ 
tudes  with  distance  to  define  and  anelastic  attenuation  model,  and  synthetic 
seismograms  for  the  derived  laterally  var)dng  model  to  estimate  the  source 
exdtation.  The  frequency  dependence  of  source  exdtation  is  considered  by 
comparing  observed  and  pre^cted  time  histories  after  being  passed  through 
narrow  bandpass  filters.  A  step  pressure  source  time  function  is  adequate  for 
explaining  the  observed  1-5  Hz  short  period  surface  waves.  The  derived 
moments  are  compared  to  other  estimates  available  in  open  literature. 

INTRODUCTION 

The  primary  focus  has  been  on  developing  techniques  for  the  estimation 
of  isotropic  moment,  M/,  from  fundamental  mode  surface-wave  data  with 
application  to  small  explosions.  Since  the  frequendes  of  interest  are  in  the  0.5 
-  5.0  Hz,  range,  path  specific  structure  and  andastic  attenuation  must  be 
defined  before  M/  can  be  estimated.  To  do  this  a  robust  technique  for  infer¬ 
ring  the  path  effects  from  a  single  station  recording  of  an  explosion  must  be 
developed  and  then  the  methodology  must  be  carefully  tested  before  being 
routinely  applied. 

We  have  made  extensive  use  of  a  USGS-GSC  refiraction  survey  in  tiie 
state  of  Maine  (Murphy  and  Luetgert,  1986;  Murphy  and  Luetgert,  1987). 
Figure  1  presents  the  profiles  analyzed.  The  profiles  are  interesting  since  they 
are  either  parallel  or  perpendicular  to  the  structural  trend  of  the  Appalachian 
Mountains  in  the  state. 

We  noted  significant  lateral  changes  in  the  moveout  of  the  surface-wave 
train.  Figure  2  shows  some  data  acquired  from  Shot  Point  6,  in  a  direction 
toward  Shot  Point  7.  The  observed  data  have  been  bandpass  filtered  between 
0.5  and  5.0  Hz.  Lateral  variation  is  seen  in  the  observed  surface-wave  train, 
center,  and  also  in  the  S-arrivals,  bottom.  While  some  of  the  S-wave  arrival 
lateral  variation  could  be  attributed  to  an  increase  of  velocity  with  depth,  it 
occurs  at  the  same  places  where  the  surface-wave  moveout  changes.  The 
arrivals  beyond  30  are  delayed,  and  cannot  be  explained  by  refracted 
arrivals.  Other  shot  also  show  a  strong  correlation  between  the  lateral  varia¬ 
tion  in  the  surface-wave  and  S-wave  arrivals. 
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Rg.  1.  Location  of  shotpoints  and  obsarvation  points  for  the  Maine  refraction 
profiles. 

1^  udng  a  sequence  of  staddng  techniques,  visual  examination, 
wavefrarm  modding  and  cross<orrdation  techniques,  a  bloddike  velocity 
modd  was  derived  (KeiUs-BOTolc  et  d,  1989;Ndet,  1990).  Table  1  shows  the 
modd  devdoped  for  this  profile.  This  table  shows  the  width  of  each  Mode, 
and  the  shear-wave  velocity  values  within  each  blodc.  The  same  Q  modd  is 
assumed  for  each  separate  vdodty  Mock.  This  model  is  shown  graphically  in 
Figures. 

ISOTROnC  MOMENT  ESTIMATES 

One  test  of  the  modd  is  to  compare  both  dte  shape  and  amplitudes  of 
(^thetk  time  hiatcMries.  The  synthetk  time  histories  are  show  at  the  top  of 
Rgure  2.  The  arrival  times  and  waveform  shapes  are  well  modded.  Another 
test  is  to  compare  peak  amplitudes  of  the  syntiketics  and  observed  data  in  the 
time  dmnain.  This  is  shown  in  Figure  4  where  the  peak  amplitudes  resulting 
frmn  different  dioices  of  bandpass  filters  are  plotted.  In  this  figure,  the  indi- 
vidud  filter  responses  are  shifted  vertically  for  clarity.  The  synthetic  data  are 
given  fay  sc^  curves  and  the  observed  data  are  given  by  the  different  sym¬ 
bols.  This  ^pire  shows  that  the  Q  modd  of  Table  1  can  explain  the  spatial 
decay  of  am]ditudes.  One  final  chedk  is  the  test  of  the  as8umpti<m  of  a  step 
source  time  function  used  in  the  sjmthetics,  by  computing  m,  for  each 
bandpass  filter.  If  the  assumption  of  the  source  time  function  shape  is 
correct,  then  the  isotropic  moment  estimates  from  the  different  bandpass 
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Fig.  2.  Seismograms  from  Shot  Point  6.  Top)  synthetic  surface-wave 
seismograms;  Center)  observed  seismograms  displayed  with  a  reduction 
velocity  of  2.5  km/s;  Bottom)  observed  shear  wave  arrivals  displayed  with  a 
reduction  velocity  of  3.5  km/s. 


Tabic  1. 

Fook  VcM^  Mn4ite  andQc  Sbuctmc  Car 
■agon  A.  Sr»  FocMm-SUc  (SQ  Obtained  fitom  Wavcfonn 
Imcntoiw  Two  Station  Fhaac  Ad^iiatBAcnt  FoOmrad  by  Picpcriion  Inversion. 
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Rg.  3.  Rrt— ntitkm  ol  the  ckrived  btenily  varying  model.  For  each  ^e 
aha^ni  tepraaafilt  the  percent  deviation  of  the  layer  velocity  fibaa  an 
ufiwriSfkted  mean  of  the  latendy  varying  vdodty. 


filtered  data  sets  should  agree.  Table  2  shows  this  comparison.  The  assump¬ 
tion  of  a  step  source  time  function  is  valid  in  the  frequency  range  of  0.5-5.0 
Hz  Ux  this  data  set.  Higher  or  lower  frequency  information  was  not  resolved 
in  the  data. 


E 

< 


> 


Fig.  4.  Comparison  of  amplitudes  of  (4»erved  (symbols)  and  synthetic  time 
histories  (curves)  after  bei^  filtered  through  different  baxulpass  filter  bands. 
The  frequencies  shown  are  cmmers  of  a  ^tterworth  bandpass  filter,  which 
has  12dWoctave  roll  off  outside  the  passband. 


We  were  aUe  to  perfrnrm  this  analysis  on  the  data  frrom  eight  shot  points, 
estimating  M/  and  aim  frnr  each  ^ot.  The  was  estimated  using  the 
density  and  compressional  wave  vdodty  at  the  source  depth  inferred  from 
the  sorfene-wave  inversion.  The  results  of  this  analysis  are  compared  to 
other  in  published  Uterature  in  TaUe  3.  The  comparison  is  also  shown  in  Fig¬ 
ures  5  and  6.  Several  features  stand  out.  First,  the  present  data  set  of  chemi¬ 
cal  cxfrioskms  indndes  only  small  yMd  events.  Se<»rul,  a  trend  for  chemical 
caqrioskms  is  not  obvious  due  to  the  differeiu»s  in  the  ejqpoiments,  e.g., 
adtether  shot  point  physical  parameters  were  known,  etc.  The  possifete 
f^tionsh^  between  the  chemical  and  nudear  exjdosion  populations  is  not 
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Table  2.  Estiaietion  of  Mi  from  filtered  seismoerams 


apparent  given  the  lack  oi  ovcriap. 


Y  (TONS) 


Fig.  5.  Con^fwriaon  of  the  inoment*to-yidd  ratio  as  a  function  of  yidkl  f(» 
eqdoakma. 

CONCXUSIONS  AND  RECOMMENDATIONS 

From  die  Mandpoint  of  wave  {Mro^mgation,  the  lateral  variations  observed 
are  larger  aome  as  large  as  10%.  Anatysis  and  modeling  of  these  data  r«{uire 
new  tedmlqiies. 

botropic  mcMnent  estimates  show  variability  f(»r  nominaDy  the  same  shot 
size.  Furtticr  refiiwment  is  not  possiUe  since  detailed  infnmaticm  on  shot 


Y  (TONS) 


Fig.  6.  Comparison  of  the  ^(o>)-to-yield  ratio  as  a  function  of  yield  for 
explosion. 

emplacement  and  shot  medium  properties  are  not  available  to  us.  On  the 
other  hand  these  results  good  since  the  multichannel  processing  gave  reliable 
estimates  of  dispersion  and  anelastic  attenuation. 

To  further  understand  the  relation  between  chemical  and  nuclear  explo¬ 
sion  yield  estimates,  it  is  necessary  to  obtain  data  in  the  10  -  100  ton  range. 
A  coordinated  field  experiment  to  acquire  such  data  from  the  shot  point  out 
to  300  km  should  be  sponsored. 

Finally,  the  data  sets  considered  indicate  that  a  step  source  time  func¬ 
tion  is  appropriate  for  the  moment  tensor.  For  large  chemical  explosions  that 
are  part  of  an  industrial  process,  the  distributed  nature  of  the  source  together 
with  the  time  delays  us^  in  blasting  would  be  expected  to  give  a  frequency 
dependent  seismic  moment  estimate  in  the  1-5  Hz  frequency  band  analyzed 
hete.  If  this  is  true,  then  the  source  spectrum  shape  wo^d  provide  a  discrim¬ 
inant  between  extended  and  point  sources  in  addition  to  the  currently  known 
ripple  firing  scallopii^  effect  on  body  wave  arrivals.  Note  that  the  surface 
wave  technique  focuses  on  the  q>ectral  shape  within  a  narrow  frequenqr 
bend  and  does  not  have  the  bwd  width  necessary  to  identify  spectral 
minima. 
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